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Q . Abstract. Supernova remnant (SNR) blast shells can reach the flow speed Vg — 0.1 c 

^ \ and shocks form at its front. Instabilities driven by shock-reflected ion beams heat the 

52 ' plasma in the foreshock, which may inject particles into diffusive acceleration. The 

1—^. ion beams can have the speed Vh ~ Vg. For Vb <C Vg the Buneman or upper- hybrid 

, ' instabilities dominate, while for Vh 3> Vg the filamentation and mixed modes grow 

^ ■ faster. Here the relevant waves for v^ « Vg are examined and how they interact 

f"^ , nonlinearly with the particles. The collision of two plasma clouds at the speed Vg is 

J-_»^ ' modelled with particle-in-cell (PIC) simulations, which convect with them magnetic 

fvj . fields oriented perpendicular to their flow velocity vector. One simulation models 

l ' \ equally dense clouds and the other one uses a density ratio of 2. Both simulations 

f~^ ' show upper-hybrid waves that are planar over large spatial intervals and that accelerate 

OO , electrons to ~ 10 keV. The symmetric collision yields only short oscillatory wave pulses, 

^7^ ' while the asymmetric collision also produces large-scale electric fields, probably through 

K^ ■ a magnetic pressure gradient. The large-scale fields destroy the electron phase space 

holes and they accelerate the ions, which facilitates the formation of a precursor shock. 
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1. Introduction 

It is well-establislied tliat a supernova remnant (SNR) can accelerate particles to energies 
~ 10^^ eV, probably through diffusive acceleration [1-3] by shocks between the SNR 
blast shell and the interstellar medium. In particular perpendicular shocks, for which 
the ambient magnetic field vector is orthogonal to the shock normal, are likely to be 
efficient particle (i.e. electron or proton) Fermi accelerators [1]. Diffusive acceleration 
across perpendicular shocks requires that the particle gyroradius exceeds the shock 
thickness. In the case of the electrons, suprathermal populations with energies in excess 
of 100 keV |:5j must be present to start with the acceleration. This is the so-called 
'injection problem'. The injection of electrons is probably achieved by waves in the 
foreshock and by the shock precursors [6-13]. 

Shocks are sources of ion beams, which originate from upstream ions that have 
been reflected by the shock or from downstream ions that leaked through the shock 
[HI [TU [15]. The spectrum of beam-driven waves depends on the beam speed vi, = \vb\ 
in the upstream plasma. The typical values of vi, are comparable to the shock speed 
Vc = |vc|, which is Vc < 10"^ c for most solar system plasmas. The corresponding Vb 
imply the growth of quasi-planar waves by Buneman instabilities, either if no magnetic 
field Bq is present or if v^ || Bq. The Buneman instability is modified if v;, ■ Bo 7^ and 
upper-hybrid waves grow if v;, ± Bq. Both waves saturate by forming electron phase 
space holes [16-19]. Upper-hybrid waves accelerate further the trapped electrons by their 
transport across the magnetic field, by which they form high-dimensional phase space 
structures [20]. Higher shock speeds of | v^ > 0.5 c exist at some SNRs and at gamma ray 
bursts [2T1[22] and ion beams moving at such speeds give rise to filamentation- and mixed 
mode instabilities [23-27]. Both modes saturate by forming current channels, which 
thermalize the plasma by their mergers. The non-linear evolution of the instability 
and the interaction of the waves with the particles depends on whether it is the 
Buneman instability, the quasi-electrostatic upper-hybrid instability or the mixed mode 
and filamentation instabilities with their electromagnetic component. The type of the 
instability thus determines the electron acceleration and injection efficiency. 

Shock speeds of Vc ~ 0.1c, which are associated with the fastest supernova-driven 
bulk flows [2n [28] . are between both extremes and it is not a priory known which 
instabilities are driven by ion beams with v^ ~ 0.1c. Wave instabilities driven by 
such beams have been examined with one-dimensional PIC simulations [12] and two- 
dimensional ones [13] with periodic boundary conditions, which limits the resolved wave 
spectrum and neglects the effects introduced by the beam front. The higher-dimensional 
wave spectrum driven by plasma beams with Vb ~ 0.1c, with a finite beam length and 
with a perpendicular magnetic field has remained practically unexplored, in spite of the 
importance of such a system for the electron injection at SNR shocks. Here we perform 
two simulation case studies that address such a system and we compare their results. 

We perform 2(1/2)D PIC simulations of colliding plasma clouds in the x — y plane, 
which evolve in time the three-dimensional particle momenta and electromagnetic field 
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vectors. SNR shocks are typically magnetized and the ratio between the electron plasma 
frequency and the electron cyclotron frequency may be a few ten to hundred. Plasma 
processes close to the front of SNR blast shells can amplify the magnetic field in the 
upstream to give a lower ratio [29l [30] . An uniform magnetic field along the z-direction 
is introduced as an initial condition. Its strength gives a frequency ratio of 5. This low 
value allows us to examine with PIC simulations the magnetic field effects |20|. Both 
colliding plasma clouds have the same density in one simulation, which gives a symmetric 
collision. A density ratio of 2 is employed in the second simulation. The spectrum of 
unstable waves is a function of the relative beam density and we can examine how robust 
the wave spectrum is against changes in it. The speed, with which both clouds collide 
along the x-direction, is set to |vc| = 0.1 c and the ion beam instabilities developing out 
of the collision have avb = f c- The electrons are heated and reach peak values of 20 keV 
in the cloud overlap region, which expands in time due to the ion motion. This value 
is in line with the results of previous PIC simulations of spatially homogeneous beam 
instabilities in the presence of a perpendicular magnetic field ^2i [13] . 

The simulations show that the dominant wave structures are quasi-planar. Two- 
dimensional electromagnetic field structures develop well behind the ion beam front. 
They fluctuate and they have a lower amplitude than the planar electrostatic fields. 
Some structuring orthogonal to the flow velocity vector is observed in an interval close 
to the ion beam front, but the electric field amplitudes remain below those of the planar 
electrostatic fields. Non-planar structures may, however, form on larger scales |31j . 
This planarity of the wave structures demonstrates that previous ID PIC simulations 
[H [H [To] should be adequate for the modelling of SNR shocks with Vc < 0.1c. 

The magnetic field is compressed in the overlap region of both plasma clouds. The 
magnetic field amplitude doubles and it is constant in the overlap region, if both clouds 
are equally dense. The collision of the asymmetric clouds amplifies the magnetic field 
close to the front of the dense ion beam to over three times the initial strength and a 
large-scale electrostatic field along Vc develops, probably through the magnetic pressure 
gradient. It accelerates ions in the flow direction and increases their speed by 20%. 

The structure of this paper is as follows. The section 2 introduces the PIC method 
and the initial conditions. The section 3 presents the results of both simulations. These 
results are discussed in the context of supernova remnant shocks in section 4. 

2. The particle-in-cell simulation and the initial conditions 

2.1. Initial plasma parameters 

We consider a system with the size L^ x Ly in the x, |/-plane. Two plasma clouds 
collide along the x-direction, each consisting of the electrons and ions. Each cloud is 
spatially homogeneous within the interval it occupies. The plasma cloud 1 is located in 
the interval —Lx/2 < x < and it moves at the speed Vb = 0.05c to increasing values 
of X. The plasma cloud 2 in the interval < a; < Lx/2 moves at the speed —Vb- The 
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relative flow velocity vector is v^ = {2vb, 0, 0). The boundary conditions are periodic in 
all directions, i.e. no further plasma enters through the boundaries. 

The plasma frequency of the electrons of cloud 1 is cJe,! = (e^ne/eome) , where 
e, nif. and Ue are the elementary charge, the electron mass and the electron number 
density, respectively. The ion mass nii = RrUe with R = 10^. The ion plasma frequency 
of cloud 1 is uji^i = uJe,iR~^^'^- The density of the plasma cloud 2 equals that of cloud 1 
in the case study 1 (run 1) and the electron and ion plasma frequencies of cloud 2 are 
^e,2 = ^e,i and Ui^2 = ^i,i, respectively. In the case study 2 (run 2) we set uje,2 = ^e,i/V^ 
and uji2 = uJi^i/\/2. The initial temperature of all plasma species in both simulations is 
T = 1.16 X 10^ K (10 eV). We use the electron skin depth A^ = c/uJe^i of the cloud 1 as 
the length scale. A spatially uniform initial magnetic field Bo = (0, 0, -Bz,o) is applied in 
both simulations with Uce/uJe,i = 0.2, where the electron gyrofrequency Uce = ^Bz^o/me- 
The cloud motion yields a convection electric field Ec = \vb x Bo|. The initial electric 
field is Eo = (0, E^, 0) for x < and Eq = (0, -E^, 0) for x > 0. The electron 
Debye length of the dense cloud is Xd,i = As/226. The electron thermal gyroradius is 
''^th,e ~ As/32. The gyroradius of the ion beams is r^ j ^ 125As for a mean speed Vb- 

The plasma clouds interpenetrate and instabilities thermalize the electrons. After 
some time, the system consists of an ion beam that interacts with the beam of incoming 
electrons. This would be a Buneman instability in the absence of a magnetic field and an 
upper-hybrid instability for the Bq we use [I2], provided that the mixed modes and the 
filamentation modes do not outgrow these electrostatic waves [27]. The fastest-growing 
electrostatic waves have a real value of the frequency in the rest frame of the electrons 
that equals the electron plasma frequency (See [12] and the cited references). This would 
be Uu = cUe,! in the reference frame of cloud 1 or Uu = oje.2 in the reference frame of cloud 
2. The wavenumber of the most unstable wave is k^ = uJuf^Vb since the speed difference 
between the ion beam and the incoming electrons is 2vb- Its wavelength A„ = 27r/fc„. 
Note that the initial conditions, in particular the spatially uniform distributions, are 
idealized. A more realistic system is discussed by Ref. [32]. However, the uniform 
beam velocities allow for a direct comparison with previous PIC simulations of spatially 
unbounded beams [12], by which similarities and differences are easily detected. 

2.2. The simulation code and resolution 

A particle-in-cell code approximates a collision-less plasma by an incompressible phase 
space fluid. The phase space fluid is represented by a large number of volume elements, 
which we call computational particles (CPs). Each CP denoted by the index i of species 
j has a charge Qj and mass ruj. It has the same charge-to-mass ratio as the particle 
species it represents, e.g. Qj/mj = —e/nie for the electrons. The absolute values of the 
charge and mass can be much larger. Each CP evolves in time under the influence of the 
collective electromagnetic fields that have been computed from all CPs. The plasma is 
updated by the code with the help of the Maxwell equations and the relativistic Lorentz 
equation of motion for each CP. The equations can be normalized (subscript A^) by 
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the substitutions E = cmeCJe, iEat/c and B = meCUe^iBjv/e for the fields. The positions 
and times are x = A^x^v and t = tj^/ujf.^i and the differential operators V = A'^Vat 
and dt = dtN/^e,i- The velocity and momentum transformations are v = cv^ and 
p = cmePAT. The masses and charges of the individual plasma species are m^ = mgrn^vj- 
and qj = eqNj and p = UeepN- The normalized equations are then 

V^ X Etv = — -^, Viv X Btv = Jtv + ^^, (1) 

OtN OtN 

^ = ^^(E^ + v^,xB^), (2) 

Vn-^n = Pn, Vn-Bn = 0. (3) 

The normalized fields Ejv and Bat can be converted into physical units, once uJe,i 
is given. The electric and magnetic fields are defined on a spatial grid, while the CPs 
follow continuous trajectories. The position of the nodes, at which the fields are defined, 
are placed in our code at equal distances. The considered system size is L^^ = 80 A^ and 
Ly = 3 Xs and it is represented by 9000 x 340 simulation cells. Each of the four plasma 
species is represented by 64 particles per cell (PPC). We use normalized units and drop 
the subscript A^. The time step A^ ?s 4.33 x 10"'^. The simulation time Tt = 225 during 
which each cloud propagates the distance 1270Aa; along the x-direction. 

3. Simulation results 

The relevant simulated field- and particle phase space distributions are examined at 
selected times. Changes in the E'^-field and in the p^-components are small due to 
the chosen system geometry and these components are not analysed. The considered 
electric E^ and Ey components are then compared to their time-animations during 
11.25 < t < 225 and in the interval y < Ly/2 . The movies 1 and 2 show E^/E^ and 
Ey/ Ef. of the run 1 and the movies 3 and 4 show E^/ E^. and Ey/ E^. of the run 2. 

3.1. Collision of equally dense clouds: run 1 

Figure [1] displays the electric E^ and Ey components at t = 56. Both plasma clouds 
overlap in the interval —2.8 < x < 2.8, in which the Er^ rapidly oscillates. The 
oscillations of Ey around Ey = have a much lower amplitude. This is in line with what 
we expect from the quasi-electrostatic upper-hybrid waves, because the electrostatic 
wave component points along the beam velocity vector. This observation is also in 
agreement with the movies 1 and 2 that initially show planar wave structures, which 
are polarized in the x-direction. However, the observed wavelength A of the oscillations 
close to the ion beam fronts at |a;| ~ 2.8 is less than the expected A„ ~ 0.6As. 

The movie 2 reveals planar wave structures in the E'^^-component. These waves are a 
consequence of the initial conditions. The abrupt change of the convective electric field 
at a; = 0, which is polarized along the y-direction, can couple to a wide k-spectrum 
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Figure 1. The electric fields E^ (black) and Ey (blue) at the position y = and time 
t — 56: The electric fields are practically anti-symmetric to x = 0. The E^ field is 
oscillatory in the interval —3 < a; < 3, in which the clouds overlap. The Ey in this 
interval oscillates with a low amplitude around Ey = 0. Wave forerunners with an 
exponentially decreasing amplitude modulus \E,j;\ are found for |a;| > 3. 
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Figure 2. The projected electron phase space distributions in run 1 at the time 
t = 56: The panel (a) and (b) show the 10-logarithm of the x, v^ and x, Vy distributions, 
integrated over < y < 0.1 and in units of a computational electron. The electrons 
are accelerated in v^ and Vy by their interaction with the upper- hybrid waves and Bq. 



of waves with k along the x-direction. These perturbations propagate in the fast 
extraordinary mode. The mean value < Ey >= in the interval \x\ < 2.8, which 
is the convective electric field, evidences that there is no net flow of plasma across Bq 
in the region in which the plasma clouds overlap. This had to be expected due to the 
equal cloud densities. Outside |x| < 3 the modulus of E^ decreases exponentially and 
the Ey component converges to the value of the convective electric field. 

Figure [2] shows the corresponding electron phase space distribution. The electrons 
have been accelerated in the interval —2.8 < x < 2.8, in which the E^ component 
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Figure 3. The B^/Ezfi component sampled at y = 0: The curves represent the time 
t = 56 (red), t = 112 (blue), t = 168 (black) and t = 225 (green). The magnetic field 
is increased from the initial value to twice that value in the cloud overlap region. The 
overlap region with the amplified Bz expands in time, due to the ion beam propagation. 



oscillates with a high amplitude. The electrons are modulated coherently in the x, Vx 
plane only close to the respective cloud fronts. The modulation does not show the 
closed elliptic orbits, which are a characteristic of the saturated waves driven by the 
Buneman instability [33] and which can initially be observed also for saturated upper- 
hybrid waves [20]. The phase speed of these waves equals the ion beam speed. The 
electrostatic potential, which traps the electrons, should thus move in the reference 
frame of the simulation box at v^ for x > or — f^ for a; < [33]. The orbits of the 
trapped electrons are, however, not centred at these speeds in Fig. [2](a). 

Instead, the centres of the phase space structures in Fig. [2]^a) move at f ~ 0. 
We attribute this and the absence of closed electron orbits to the magnetic field that 
accelerates the electrons in the fy-direction in Fig. Ej^b). This acceleration imposes a 
drag on the wave that is responsible for the reduced phase speed [31], explaining also 
that the wavelengths of the oscillations of E^ in Fig. [1] are shorter than A„. Note, 
however, that the electrons rapidly detrap and that they do thus not form the phase 
space cones [20] that evidence electron surfing acceleration. The phase space structures 
in the x, Vx plane in the interval —1.5 < x < 1.5 show twice the length of those at the 
cloud fronts. They have been formed by the coalescence instability [16]. Figure [2](b) 
also reveals that cusps in the Vy speed of the electrons have developed at the leading 
edges of the clouds at x ~ ±2.8, which must be linked to a B-^ that varies along x. 

The magnetic B^ component is indeed not spatially uniform. Figure [3] demonstrates 
that the value of the background magnetic field is amplified to B^/Bz^ = 2 in the 
cloud overlap region. Snapshots at four simulation times evidence that the interval 
with Bz ^ 2 Bzfl expands at a constant rate. This magnetic expansion is driven by 
the ion beams that do not change their speed much during the simulation time. The 
magnetic field amplitude will increase beyond 25^0 when a shock starts to form, as PIC 
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Figure 4. The electric fields in simulation 1 at t = 225: Panel (a) shows the 
Ex component, which is due to the upper- hybrid wave with its k || v;,. It shows 
planar structures. Panel (b) shows Ey. The oblique mixed modes and charge density 
modulations by electron bunches contribute to this electric field component. 



simulations of perpendicular shocks demonstrate [8], [9l [TO] . 

The movie 1 shows that the electric field is predominantly polarized along the x- 
direction and that the field structures are planar before t ^ 56. After this time the 
Ey-S.e\d component starts to be structured along the y-direction, which is usually a 
characteristic of oblique modes and of filamentation modes. It can, however, also be 
caused by bunches of accelerated electrons that yield charge density fiuctuations. Figure 
m shows the E^ and Ey components at t = 225. Both fields have amplitudes exceeding 
Ec- The E^ is strongest and shows planar wave fronts, in particular close to the leading 
edge of the cloud at x ~ 11. Filamentary structures in Ey are observed in the overlap 
region of both clouds, but their amplitude is lower. The strongest oscillations of Ey 
reach 2 Ec at the leading edge. The E^ oscillations are 2.5 times stronger than those 
of Ey. The movie 2 shows that the structuring of Ey at the ion beam front along the 
y-direction changes only slowly in time. The oscillations of Ey at 10 < x < 12 in Fig. 
m^b) are thus caused by an instability between the ion beam and the incoming plasma 
[27] . The structures for x < 8 in Fig. \Mjo) fiuctuate according to the movie 2 and they 
are probably caused by fast electron bunches that move in the x-y plane. 

Figure [5] shows the final phase space distributions of the electrons and ions. We 

1 /9 

compute the electron phase space density as a function of x and v± = {vl + Vy) and 
the ion phase space density as a function of x and v^- The distributions are integrated 
over the interval < y < 0.1. We find that the electrons are accelerated to about 5vf, 
(16 keV), which is 2.5 times the streaming speed of the ions relative to the incoming 
electrons. The momentum modulation of the ion beam is of the order of 1% and rises to 
about 10% at the cloud front. The electrons in the overlap region show a spatially almost 
uniform velocity distribution. The ions, on the other hand, have not yet thermalized. 
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Figure 5. The particle distributions in nm 1 integrated over < y < 0.1 at i = 225: 
Panel (a) displays the 10- logarithm of the electron distribution, as a function of v±. 
The electrons are rapidly heated by the ion beam front. Panel (b) shows the 10- 
logarithmic ion beam distribution, which shows phase space holes at its front. Both 
panels show the phase space density in units of a computational particle. 
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Figure 6. The electric fields E^ (black) and Ey (blue) at the position y — and time 
t ~ 56: The E^ field is oscillatory in the interval — 3 < x < 0. A single wave structure 



occurs in the Ex-Hcld for < a; < 5. Its oscillation period far exceeds A^, 
this interval oscillates with a low amplitude around Ey = 0. 



The E„ in 



3.2. Collision of clouds with non-equal densities: run 2 

Figure [6] displays the electric E^ and Ey components at t = 56. The E'^.-component 
behaves qualitatively different in the interval x < 0, in which the tenuous ion beam 
interacts with the dense plasma of cloud 1 and in the interval x > 0, where the dense 
ion beam interacts with the tenuous plasma of cloud 2. The oscillations for x < have 
a A ~ A„/2, as in run 1. Only weak oscillations with A ~ A„/2 are observed at x ~ 3, 
which are superposed onto a powerful wave structure with A ^ A^. This electric field 



Simulations of ion beam instabilities 



10 




X-Position 



Figure 7. The Bz/Bzfi sampled at y = 0: The curves represent the time t — 56 (red), 
t = 112 (blue), t = 168 (black) and t = 225 (green). The B^ is amplified in the cloud 
overlap layer. The maximum of B^ is located close to the front of the dense ion beam. 
It propagates with the beam and the maximum of B^ increases with time. 



forms a ramp for the ions. The ions of cloud 2, which move to decreasing values of x, will 
be slowed down in the interval 1 < x < 5 and the ions of cloud 1 moving to increasing x 
are decelerated in < x < 1. The movie 3 demonstrates that the large-scale structure 
is long-lived and that it does not develop a significant structuring in the y-direction. 
The movie 4 shows that this large-scale structure has no Ey component. 

The large-scale electric field structure is correlated with a change in the strength 
of B^. The Bz has been compressed during the time interval up to t = 56. Figure [7] 
demonstrates this for y = sampled at four times. The Bz a.t t = 56 peaks at x ~ 1.2, 
which coincides with the position where the E^ of the large-scale electric field structure 
in Fig. [6] changes its sign. A comparison of E^ with 5^ at t = 56 suggests that Er^ 
changes with dBz/dx. This can be interpreted as E^ oc Bz{dBz/dx), since B^ > 
everywhere. The exact proportionality cannot be demonstrated, because short waves 
are superposed on the large-scale structure. A proportionality of E^ oc Bz{dBz/dx) 
would imply that the large-scale E^ is driven by the magnetic pressure gradient. Figure 
[7] also reveals that the amplification of Bz is lower for x < in comparison to run 1. 

Figure [H] displays the electron distribution in run 2 at t = 56. The x, Vx distribution 
shows well-defined electron phase space holes that involve the bulk of the electrons, but 
only at the leading edge of the tenuous cloud at x ~ —2. The periodic structures in x, Vx 
are apparently stable only in the interval, in which the Vy speed of the bulk electrons 
and thus Bz does not significantly change. The rapid change in the mean Vy component 
at X ~ 1.5 destroys the electron phase space holes and the bulk electrons are heated. 
Some oscillatory electron orbits are visible for x > but the number density of the 
electrons following these paths is low. The lack of periodic electron structures with a 
significant charge density explains, why no strong oscillations in E^ on a scale ~ A„ are 
visible in this interval in Fig. [6l Like in Fig. [2t^a) the electron phase space holes in 
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Figure 8. The electron phase space distributions in run 2 at i = 56: Panel (a) 
and (b) show the 10-logarithm of the x,Vx and x,Vy distributions, integrated over 
< y < 0.1 and in units of a computational electron of the dense cloud. The electrons 
are accelerated by their interaction with the upper- hybrid wave and with Bq. 



Fig. [8](a) are tied to an electrostatic wave with a phase speed, which is lower than that 
driven by the Buneman- or the upper-hybrid wave instability. 

Figure [9] shows the E^ and Ey fields in the half-plane x < at t = 225. The 
amplitude of E^ is about twice that of Ey and it shows small scale, quasi-planar 
oscillations. These oscillations are characteristic for electron phase space holes that 
are extended along the y-direction in form of tubes. The Ey field shows fiuctuations 
along both directions. As we move away from the beam front at x ~ —11 to increasing 
values of x, the characteristic size of the fiuctuations in Ey increase and the structures 
in Ex become less planar. The E^- and Ey amplitudes are lower than their counterparts 
in Fig. m but their spatial distribution is similar. 

Figure fTOl shows the E^ and Ey in the half-plane x > 0. The small-scale oscillations 
of Ex are comparable to those oi Ey. A planar negative electric field structure is found 
at X = 4.5 and the oscillations close to the cloud front at x = 11 are superposed on a 
structure with a positive E^. The Ey fields are strong only in the interval 5 < x < 10. 
The large-scale fields exceed the oscillations in E^ and Ey in their amplitude and should 
be the strongest accelerators in the plasma system, due to their spatial coherence over 
more than one A^. The field structure is also stable, as the movie 3 evidences. This 
structure is again related to the steepest gradients of Bz at t = 225 in the Fig. [71 

Figure [11] shows the final electron and ion distributions in run 2, which are both 
asymmetric relative to x = 0. The electron acceleration in x < is somewhat stronger. 
The electron peak speed is comparable to that in run 1 in Fig. [51 The accelerating 
electron beams in the intervals < x < 5 and 12 < x < 15 are related to the magnetic 
field gradients in Fig. [71 and involve the Vy component. The ion beam going to positive 
X in Fig. [TTlfb) confirms what we anticipated in the discussion of Fig. [TOl The ions 
slow down in the negative electric field at x ~ 5 and accelerate at x ~ 11. The ions are 



Simulations of ion beam instabilities 



12 




Figure 9. The electric fields in simulation 2 for x < 0: Panel (a) shows the E^ 
component and panel (b) the Ey component at t = 225. The planar structures of E^ 
at cc < —10 are due to the upper-hybrid waves. These waves damp out for increasing 
X > — 10 and they loose their planarity. The typical size of the structures of -Ej^ increase 
for increasing a; > — 10 and their amplitude decreases. 
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Figure 10. The electric fields in simulation 2 for x > at t = 225: Panel (a) shows 
the Ex component and panel (b) the Ey component. Planar wave structures in E^ 
are superposed on a large-scale electric field modulation. The Ey component shows 
structuring along y at x « 9, which decays for decreasing values of x. 



accelerated by more than 20%, which exceeds the ion acceleration in run 1. Figure [TlT c) 
shows that the tenuous ions going to negative x are almost un-accelerated. The ions 
spread out at the cloud front due to their electrostatic interaction with the electrons. 
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Figure 11. The particle distributions in run 2 integrated over < j/ < 0.1 at 
t = 225: Panel (a) displays the lO-logarithm of the electron distribution, as a 



function of v 



{vl 
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Panels (b) and (c) show the 10-logarithmic phase 



space distributions of the respective beams. Both panels show the phase space density 
in units of a computational particle of the dense plasma cloud. The electrons have 
been thermalized for —8 < a; < and 5 < a; < 10 and a beam exists at a; w 4. Only 
the ions of the dense beam in (b) are significantly accelerated. 



4. Discussion 

In this work we have examined with particle-in-cell (PIC) simulations the collision of 
two magnetized plasma clouds. Plasma instabilities form in the layer, in which both 
clouds overlap. After a short transient phase the plasma thermalization is accomplished 
by an ion beam moving through a cool electron plasma. The cool electron plasma is 
from that cloud that has not yet interacted with the second cloud. The electrons of the 
second cloud have been heated. They are confined by the perpendicular magnetic field 
to the layer, in which both clouds overlap. This layer expands. One ion beam-driven 
instability develops at each of the two fronts of the overlap layer. 

The plasma clouds have been equally dense in one simulation, which is typically 
assumed by PIC simulations of kinetic shocks that employ the piston method [H [9|[T0]. 
The piston method models one plasma cloud, which it reflects at one simulation 
boundary onto itself. Two clouds are modelled here and we get two instabilities with 
equally dense ion beams, one at each front of the expanding overlap layer. The ion 
beam density is comparable to that of the electron beam, which is the equivalent of 
the Buneman instability in magnetized plasma. The density ratio has been set to 2 in 
the second simulation, which gives one instability driven by a dense ion beam and one 
driven by a tenuous ion beam. We could thus assess the impact of the ion beam density 
on the plasma dynamics for three cases using two simulations. 

The simulation parameters are probably relevant for the foreshock regions of the 
fastest perpendicular shocks at supernova remnants (SNRs). The normalization of the 
Maxwell-Lorentz set of equations makes the results independent of the value of the 
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plasma frequency in physical units. The strength of the magnetic field resulted in a 
ratio of 5 between the electron plasma frequency of the dense cloud and the electron 
cyclotron frequency. The equivalent value close to SNR shocks is thought to be several 
hundred. However, a substantial decrease of this ratio to about 10 might occur at some 
SNR shocks [29] . Shock-accelerated particles can also amplify the magnetic field locally 
[30] . We have selected a ratio of 5, because it allows us to resolve magnetic field effects 
such as the electron surfing acceleration |20j . 

Some upstream ions are reflected by the shock [6-10] and they gyrate in the 
foreshock, which is a scenario analogous to that we have modelled here. The ion beam 
speed can reach twice the shock speed in the upstream frame of reference, if the reflection 
is specular. Our ion beam speed in the reference frame of the electron beam is 0.1 c, 
which can be achieved with realistic SNR shock speeds [21], [28] . It is between the low 
speeds, for which electrostatic instabilities dominate, and the fast speeds that favor 
mixed- or filamentation modes [23l [27] . Our simulations here have shown that the 
wave structures are quasi-planar and the electrostatic instabilities thus dominate. The 
emergence of filamentary structures indicates, however, that this speed is probably the 
most extreme case in this respect, at least for our plasma parameters. 

We have introduced a second spatial dimension and an ion beam with a finite length. 
This has left unchanged the findings of previous one-dimensional PIC simulations of 
spatially uniform beam distributions with respect to the acceleration of electrons. The 
peak speed of the electrons is still limited to a few times the ion beam speed in the 
reference frame of the electrons. This peak energy has furthermore been independent of 
the ion beam density. The peak energy has been 10-20 keV for the ion beam speed 0.1 c, 
which falls short of the 100 keV that are required for the electron injection into diffusive 
acceleration [5|. Higher electron energies may be reached when the ions thermalize. 

If both plasma clouds have the same density, the ion acceleration is identical close 
to the fronts of both beams and it amounts to maximally Vb/10. The simulations 
have shown that the magnetic field is convected with the plasma slabs, as expected 
from the frozen-in theorem. The magnetic field strength doubles where both clouds 
overlap. We could observe well-developed electron phase space holes. The speed of 
their electrostatic potentials has been low, presumably due to the magnetic drag [34] 
excerted on the trapped electrons. The electron surfing acceleration [20j by the transport 
of the trapped electrons across the magnetic field has been weak. 

The ions at the front of the dense beam have been accelerated by Vb/5 in the 
simulation with the unequal cloud densities, while those at the front of the tenuous 
beam have only been scattered. This has originated in the asymmetric magnetic field 
compression. The magnetic field strength has increased by only 50% at the front of 
the tenuous beam while it has reached over 3 times its initial value at the front of the 
dense beam. The magnetic field strength has also varied more rapidly. This variation 
was correlated with a large-scale electric field, which may have been driven by the 
magnetic pressure gradient. It could strongly interact with the ions. This electric field 
has destroyed the electron phase space holes in the bulk plasma, probably because it 
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has removed the stable equihbrium position of the potential. 
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